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ABSTRACT 



We present an improved model for an asteroseismic diag- 
nostic contained in the frequency spacing of low-degree 
acoustic modes. By modelling in a realistic manner re- 
gions of rapid variation of dynamically relevant quanti- 
ties, which we call acoustic glitches, we can derive sig- 
natures of the gross properties of those glitches. In partic- 
ular, we are interested in measuring properties that are re- 
lated to the helium ionization zones and to the rapid vari- 
ation in the background state associated with the lower 
boundary of the convective envelope. The formula for 
the seismic diagnostic is tested against a sequence of the- 
oretical models of the Sun, and is compared with seismic 
diagnostics published previously by Monteiro & Thomp- 
son (1998, 2005) and by Basuetal. (2004). 

Key words: solar interior; solar oscillations; stellar oscil- 
lations. 



1. INTRODUCTION 



Abrupt variation in the stratification of a star (relative 
to the scale of the inverse radial wavenumber of a seis- 
mic mode of oscillation), such as that resulting from 
the (smooth, albeit acoustically relatively abrupt) depres- 
sion in the first adiabatic exponent 71 = (dlnp/d\n p) s 
caused by the ionization of helium, where p, p and s are 
pressure, density and specific entropy, or from the sharp 
transition from radiative to convective heat transport at 
the base of the convection zone, induces small-amplitude 
oscillatory components (with respect to frequency) in the 
spacing of the cyclic eigenfrequencies v n< i of seismic os- 
cillation. One might hope that the variation of the sound 
speed induced by helium ionization might enable one to 
determine from the low-degree eigenfrequencies a mea- 
sure that is directly related to, perhaps even almost pro- 
portional to, the helium abundance, with little contamina- 
tion from other properties of the structure of the star. 

A convenient and easily evaluated measure of the oscil- 
latory component is the second multiplet-frequency dif- 
ference with respect to order n amongst modes of like 
degree I: 



(1) 



(Gough, 1990). Any localized region of rapid variation 
of either the sound speed or the density scale height, or a 
spatial derivative of them, which here we call an acous- 
tic glitch, induces an oscillatory component in /S.2V with 
a 'cyclic frequency' approximately equal to twice the 
acoustic depth 
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of the glitch, and with an amplitude which depends on the 
amplitude of the glitch and which decays with v once the 
inverse radial wavenumber of the mode becomes compa- 
rable with or greater than the radial extent of the glitch. 

We report here on a model formula for the oscillatory 
contribution to the second frequency difference that is di- 
rectly related to the gross properties of the glitch, and 
compare how well it fits actual eigenfrequency differ- 
ences with corresponding fits of previously published for- 
mulae. 



2. THE FITTING FUNCTION 



For a nonrotating star with vanishing pressure gradient 
at the surface, the angular eigenfrequencies ui of adia- 
batic oscillation obey a variational principle (e.g. Chan- 
drasekhar 1963) which can be used to estimate the con- 
tribution Suj to u> from the acoustic glitch. In estimating 
the contribution from helium ionization we retain only 
the dominant glitch term (571/71 at fixed pressure and 
temperature), and we represent £71 /71 by a pair of (neg- 
ative) Gaussian functions of r with widths Ai and An, 
whose integrals are Ti and Tn, and which are centred 
about the acoustic depths t\ and t\\ of the first and sec- 
ond ionization zones of helium beneath the seismic sur- 
face r = R of the star. After converting to cyclic fre- 
quency v = u)/2%, the oscillatory component is given 
asymptotically by 



x ( M /3K I ~ 1 e- 87r2 ^ K i A " !y2 cos 2ipi 

+ k^c-^ 2 ^ 2 ^ 2 cos 2Vn) (3) 

(Houdek & Gough, 2007). We evaluate the phases ipj = 
ip(fi) and ipn = ip(fu), where lot — lot + en, by rep- 
resenting the envelope by a plane-parallel polytrope of 
index m =3.5 and adding a phase constant en to lot to 
account for the deviation of the actual envelope from the 
polytrope: 



iI>{t) = loth — (m + 1) cos 1 
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• (4) 



In equation (3), Ki = k(t{) etc, with k(t) = [1 — (m + 
1) 2 /47t 2 ^ 2 t 2 ] 1 / 2 ; also vq is a characteristic value of the 
first frequency difference Aiv whose value is half the in- 
verse total acoustic radius T of the star: v = lo /2n = 
1/2T. We have found that the ratios /? = I^An/TnAi, 
V = ti/tu and \i = Ai/An hardly vary amongst stel- 
lar models whose masses and radii vary by factors of at 
least five. We have accordingly adopted constant values, 
namely (3 = 0.45, r? = 0.7 and [i = 0.9 for solar-like 
stars. 

To obtain a complete description of A 2 v we must add an 
oscillatory contribution from the near discontinuity in the 
density scale height at the base r c of the convection zone: 



6 C:OSC v ~ A c vlv 2 (1 + 1/16tt 2 t V) 
x {cos[2i/; c + tan _1 (4-7rT ^)]} 
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with 
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(5) 
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where ip c — ?/>(t c + w _1 e c ) is the polytropic phase at 
the base of the convection zone, the phase constant e c 
being allowed to differ from en to account for the fact 
that the region of the convection zone beneath the helium 
ionization zones is not an exact polytrope; the quantity 
t is a characteristic scale over which the sound speed 
in the radiative interior relaxes to a smooth extrapolation 
from the convection zone (and which we also take to be 
constant, 80 s, for all the stars we consider). 

From these expressions it is straightforward to evaluate 
the second difference A2V, to which must be added a 
smooth term which we represent by a third-degree poly- 
nomial in z/ _1 : 



12, sm 



(7) 
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The eleven parameters An, A c , An, m, r c , en, e c and at 
are adjusted to fit by least squares the theoretical curve to 



the second frequency differences of the actual eigenfre- 
quencies of the modes. 

We compare the current seismic diagnostic (3)-(7) with 
diagnostics published previously by Monteiro & Thomp- 
son (1998, 2005) and by Basu et al. (2004). 



2.1. The seismic diagnostic by Monteiro & Thomp- 
son 

This diagnostic is based on a single-triangular approx- 
imation to £71/71 due to He ionization. The resulting 
formula for the frequency perturbation is 



sin 2 (27w/3) 



vj3 



cos(47tz/th + 2en) 
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cos(47TfT c + 2e c ) , (8) 



in which (3 is the half-width and 61 is proportional to 
the height of the triangle, from which second differences 
can be computed. We set to zero the fitting parameter 
C2, a nonzero value of which, in the original formulation, 
would account for an putative adiabatic extension to the 
convection zone produced by overshooting (which is not 
present in our solar models). By so doing, the seismic 
diagnostic (8) has the same number, 1 1, of fitting param- 
eters as the seismic diagnostic (3)-(7) when adopting the 
third-degree polynomial (7) to represent the smooth con- 
tribution. 



2.2. The seismic diagnostic by Basu et al. 

In this formulation the seismic diagnostic is that used by 
Basu (1997), which was invented by Basu, Antia and 
Narasimha (1994) for studying the base of the convec- 
tion zone. Ballot et al. (2004) and Piau et al. (2005) 
have adopted essentially the same formula. It is not di- 
rectly linked to the form of the ionization-induced acous- 
tic glitch, but instead was motivated by the oscillatory 
signature of a discontinuity in sound speed, whose am- 
plitude was modified to make it decay with frequency ap- 
proximately in the manner of the second differences of 
the eigenfrequencies. The functional form is given by 

A 2 S osc iy ~ ^61 + cos(47wth + 2e n ) 

+ (ci + cos(4twt c + 2e c ) . (9) 

Here we set ci = 0, to remove the effect of the erro- 
neous sound-speed discontinuity at the base of the con- 
vection zone and to render the functional form of the con- 
tribution at r = t c (nearly) consistent with those in the 
other two signatures (5) and (8); and we adopt the third- 
degree polynomial (7) to represent the smooth contribu- 
tion, thereby using the same number (11) of fitting pa- 
rameters as the seismic diagnostics (3)-(7) and (8), (7). 
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Figure 1. Denotation of the nine calibrated solar models 
which are used for testing the formulation of the second 
frequency differences. The ' central model' is model 0; the 
set of the four models ( 1 )-(4) have a constant value of Z 
but varying age; the second set of the models, (5)-(8), 
have constant age but varying Z. 



3. TEST MODELS 

Two sets of calibrated evolutionary models of the Sun, 
computed by Gough and Novotny (1990), were used. 
One set of models has a constant value for the heavy- 
element abundance Z but with age varying in uniform 
steps of approximately 0.051 in \nt, and the other set 
has constant age but varying heavy-element abundance in 
uniform steps of approximately 0.016 in In Z (see Fig. 1). 



4. RESULTS 

All three diagnostics were fitted by the same numerical 
algorithm. The lower panels of Figs 2-4 show the re- 
sults for the three seismic diagnostics (3)-(7), (8) and 
(7), and (9) and (7), fitted by least-squares to the adiabat- 
ically computed eigenfrequencies of the nine test mod- 
els. Plotted are the frequency-averaged amplitudes, (Au) 
and (A c ), of the two oscillatory seismic contributions as- 
sociated with Hen ionization and with the base of the 
convection zone, obtained by averaging the appropriate 
calibrated formula uniformly over the frequency range 
of the data plotted in the upper panels. Also plotted are 
the values for the acoustic depths Tu and r c . Depending 
on the seismic diagnostic, the values of either An, $ or 
the phase constant en are shown, and also the standard 
measure \ 2 of the goodness of the fit between A2f of 
the modelled frequencies and the calibrated seismic di- 
agnostic. In all three cases the frequency-averaged am- 
plitude (Au) increases with the helium abundance Y of 
the test models; moreover, the amplitudes Tn and b\ in- 
ferred from diagnostics (3)-(7) and (8), (7) respectively 
also increase with Y (Tu varying almost linearly in the 
Z-varying sequence, b\ varying almost linearly in the 
i-varying sequence), the latter result being apparently 



at variance with the findings of Monteiro and Thomp- 
son (2005). The seismic diagnostic (3)-(7) fits A2^ of 
the model eigenfrequencies the best, as is evinced by its 
smallest values of \ 2 ■ 

The quality of the fits is also illustrated in the upper pan- 
els of Figs 2-4, which show the second frequency differ- 
ences over the whole frequency range considered for the 
central model 0. It is also clear from comparing these fig- 
ures that the seismic diagnostic (3)-(7) reproduces A2V 
of the model frequencies the most faithfully. 

It is interesting to note that in determining the parameters 
in the formulae the fitting procedure gives greater weight 
to the low-frequency end of the ionization-induced com- 
ponent A2, 7 , where the amplitude is the greatest, and 
greater weight to the high-frequency end of the contri- 
bution A2, c from the base of the convection zone, where 
misfits cannot be offset by adjustments to A2. 7 ; there- 
fore, the diagnostics (8), (7) and (9), (7) fail principally 
in the mid-range. Some large-scale discrepancy can be 
seen also at high frequency for formula (9), (7) because 
A2 )7 decays with frequency too slowly. We acknowledge 
that these discrepancies might be either smaller or hardly 
greater than the random errors in imminent asteroseis- 
mic data; nevertheless, fitting an erroneous formula can 
lead to the inferences from that fitting being biased, and 
should therefore be avoided. 
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Figure 2. Top: The symbols are second differences A2V, defined by equation (1), of low-degree frequencies obtained 
from adiabatic pulsation calculations of the central model 0. The curve is the seismic diagnostic (3)-(7) whose eleven 
parameters have been adjusted to fit the data by least squares. Bottom: Model properties determined from the seismic 
diagnostic (3)-(7) for the nine test models. Models with varying age are connected by solid lines, those with varying Z 
by dashed lines. Some model identifiers are written in the upper left panel. The smallest (central) frequency value used in 
the least-square fitting is v ~ 1322 pHz, the largest is v ~ 4058 fiHz. The lower right panel is the standard measure \ 2 
for the goodness of the least-squares fit. 
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Figure 3. Top: The symbols are second differences A2V, defined by equation (1), of low-degree frequencies obtained 
from adiabatic pulsation calculations of the central model 0. The curve is the seismic diagnostic (8),(7) whose eleven 
parameters have been adjusted to fit the data by least squares. Bottom: Model properties determined from the seismic 
diagnostic (8), (7) for the nine test models. The smallest (central) frequency value used in the least-square fitting is 
v ~ 1322 fjbHz, the largest is v ~ 4058 jiHz. The lower right panel is the standard measure x 2 for the goodness of the 
least-squares fit. 
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Figure 4. Top: The symbols are second differences A2V, defined by equation (1), of low-degree frequencies obtained 
from adiabatic pulsation calculations of the central model 0. The curve is the seismic diagnostic (9),(7) whose eleven 
parameters have been adjusted to fit the data by least squares. Bottom: Model properties determined from the seismic 
diagnostic (9), (7) for the nine test models. The smallest (central) frequency value used in the least-square fitting is 
v ~ 1322 \iHz, the largest is v ~ 4058 fiHz. The lower right panel is the standard measure \ 2 for the goodness of the 
least-squares fit. 



